d = float(input('砂桩直径(m): '))
e_max = float(input('地基最大孔隙比: '))
e_min = float(input('地基最小孔隙比: '))
e0 = float(input('未加固时孔隙比: '))
d_r1 = float(input('加固后土的响度密实度: '))

e1 = e_max-d_r1*(e_max-e_min)
shape = int(input('砂井布置方式: 等边三角形(1), 正方形(2)'))
if shape == 1:
    a = 0.95
elif shape == 2:
    a = 0.89
else:
    print('请输入数字：等边三角形(1)或正方形(2):')
l1 = a*d*((1+e0)/(e0-e1))**0.5

print('正三角形布置时砂桩的中心距为 %f m' %l1)
